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Abstract 

We simulate 4d simplicial gravity for three topologies S 4 , S 3 x S 1 and S 1 x S 1 x S 1 x S 1 
and show that the free energy for these three fixed topology ensembles is the same in the 
thermodynamic limit N± — > oo. We show, that the next-to-leading order corrections, at 
least away from the critical point, can be described by kinematic sources. 



1 Introduction 



It is not a priori clear, whether in a path-integral formulation of quantum gravity the sum 
over metrics should also run over topologies. In a theory containing topology fluctuations 
only those topologies contribute to the sum in the thermodynamic limit, which maximize 
the extensive part of the free energy. Other contributions are exponentially suppressed. 
If one believes that all topology excitations should be present in the continuum theory, 
the bulk volume contribution to the free energy should be independent of topology. An 
explicit example of such a theory is provided by two dimensional quantum gravity. Here 
the Einstein-Hilbert action is of purely topological nature. It can be shown explicitly 
that the coefficent /x c of the leading (extensive) part of the entropy does not depend on 
the topology. The coefficient of the logarithmic correction is 7— 3, where 7 is the so-called 
surface susceptibility exponent. The exponent 7, does, however, depend linearly on the 
genus of the surface. The linear dependence leads to the double scaling limit at which 
one can reduce the number of non-perturbative modes of the theory with fluctuating 
topology to the solutions of the Painleve II equation. In four dimensions the situation 
is more complicated. At present, no classification of topologies is known. Therefore 
the topological part in the action is unknown. Simplicial quantum gravity allows us 
to sum over geometries with fixed topology. We investigate numerically three different 
topologies and show that in these cases, up to the leading order, the free energy does 
not depend on the topology. We observe a topology dependence in the next to leading 
order. We study those volume corrections and analyze their sources at the kinematic 
bounds. Some of our results have already been presented at Lattice 96 0. 



2 Definitions 

The partition function of simplicial quantum gravity in the grand canonical ensemble is 



where the first summation is over all 4<i-simplicial manifolds T with fixed topology 
||. The parameter K4 is proportional to the cosmological constant and k<i is a linear 
combination of the inverse of the Newton constant and the cosmological constant in the 
naive continuum limit. The prefactor 1/C(T) is a remnant of the invariance group and 
divides out the internal symmetry factor of the triangulation. The free energy F in the 
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canonical ensemble is denned as: 

e F( K2 ,N 4 ) = _J_ e ^N 2 (T) ^ 2 ) 

C(T) 

The sum runs over fixed topology Ad simplicial manifolds with a fixed number N4 of 
4-simplices. In the large volume limit N± — > oo the free energy is assumed to have the 
form : 

F(K 2 ,N i ) = NJ(K 2 ) + 5(K 2 ,N 4 ), (3) 
where the function 5 is a finite size correction, ie for any k 2 

lim = 0. (4) 

To recover some basic properties of the theory it is convenient to study the derivatives 
of the free energy. The action density 

is normalized such that it becomes an intensive quantity in the thermodynamic limit 
N± oo. It is related to the average bare Regge curvature R av by r = R av + a where 



a 
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^■arccos(l/4) »s 2.09784. The average is taken in the canonical ensemble 



The derivative with respect to N4 



dF 

k 4 (k 2 ,N 4 ) = —. (6) 

will be called the critical value of the parameter K4. Our definition differs slightly from 
the one proposed in [||]. This quantity is a measure of finite size dependence of the free 
energy, ie how fast it approaches the thermodynamic limit £4 — > K4(k 2 ) for N4 — > 00. In 
the thermodynamic limit N4 —> 00 the value £4(^2) defines a critical line of the model 
corresponding to the radius of convergence of the series (0). 

Taking the second derivative, we see that k 4 is related to the action density r in the 
following way : 

•9^4 AT dr 

a^ = r + N 'lm (7) 

It is important to note that in the large iV 4 limit the second term on the right hand side 
of (0) goes to zero, and the critical parameter £4 becomes an integral of r. Thus, if r 
is independent of topology in the thermodynamic limit, so is £4, unless the integration 
constant depends on the topology. To fix the integration constant, one just has to 
measure £4 for one particular value of k 2 . 
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3 Methods 



The average action density fl5|) can easily be measured in canonical simulations. The 
quantity k 4 ([]) requires non- conventional methods. Some of them, like those based on 
the sum rules [|J or an analysis of the baby universe distributions || , permit to directly 
extract the next-to-leading volume corrections. These methods are unfortunately limited 
to the elongated phase. Following ||, we adopt here a more general method based on 
multi canonical simulations which works equally well in the entire range of the coupling 
k 2 . To learn about how the free energy depends on N 4 one lets the volume fluctuate in 
the external potential U(N 4 ). Measuring the resulting N 4 distribution and combining 
it with the known form of U one gets the dependence of the free energy on N 4 . The 
freedom one has in choosing of the potential U should be used to minimize the error for 
the quantity one wants to measure. In this particular case one wants to measure 8F/dN 4 
for the function F, which is expected to smoothly approach a function A^ 4 /(k 2 ) linear in 
iV4 for large iV 4 as given in eq. (||), (f|). A Gaussian term controlling fluctuations around 
a fixed volume V 4 is well suited for this problem || : 

U = -k 4 N 4 + ^(N 4 -V 4 ) 2 (8) 

but other terms can be used as well 0. The multicanonical partition function for this 
potential reads : 



Z = Y^ e K2N2+u{Ni) = e FiK2 ' N4) ~ K4N/i+ 2 {Ni ~ V4)2 (9) 

T N A 

Expanding F(k 2 , N 4 ) around N 4 = V 4 one sees that the N 4 distribution becomes Gaus- 
sian : 

P(x = N 4 - V 4 ) ~ exp(-^(x - x ) 2 + ...) (10) 

when k 4 is tuned to be close to the derivative dF/dV 4 and 7 is much larger than the 
second derivative d 2 F/dV 4 2 . This means that the range of the distribution is much 
smaller than the typical scale for changes in F. The parameters V and Xq are related to 
F and U : 

r = 7 - d 2 F/dVi, x = (dF/dV 4 - k 4 )/T. (11) 

Both these quantities can be measured in the simulations, namely T from the width 
of the distribution and Xq from the shift of the maximum from N 4 = V 4 . From the 
numerical results one gets an estimator for the critical coupling at V 4 : 

k 4 = 0F/dV 4 = Tx + k 4 . (12) 
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In this formula, k 4 is the coupling used in the simulation, from which one extracted Xq 
and T. The estimator can be improved by minimizing xq. This can be done recursively 
by setting k 4 = k 4 in the potential . 

After performing measurements one can justify the validity of the assumptions used to 
write the Gaussian approximation fllCf ). The value for 7 should not be too large, because 
suppressing volume-fluctuations can spoil the mobility of the algorithm. On the other 
hand, it should not be too small either, because the average number of sweeps between 
two configurations with canonical volume V4 grows with decreasing 7. In particular 
we checked that for 7 = 0.0001, T computed from the width of the distribution, T = 
l/y / o" 2 (A^4), i.e. AN 4 = 100, was equal to 7 within the error bars, meaning that 7 ^> 
d 2 F/d 2 V4. Therefore the free energy changes very slowly in the range of the distribution 
width as needed for the approximation fllPl) . 

In this investigation we have performed simulations with V4 = 4000, 8000, 16000, 
32000, 64000. 

Because the local moves used in the Monte Carlo simulation preserve the topology t 
of the manifold, t is given by the topology of the starting configuration. In our simulation 
we choose to simulate a spherical S A topology and two different tori, namely S 3 x S 1 
and ^xS'x^x S 1 . 

The starting configuration for the sphere is, as usual, the 4<i-boundary of a 5-simplex. 
The manifold S 3 x S l can be produced from a spherical manifold by taking away two 
separate 4-simplices and gluing together the boundaries 61,62, which are created when 
removing the two 4-simplices. However, two vertices in a 4d simplicial manifold may 
have at most one common link. To avoid creation of a double connection in the gluing 
procedure, the distance between the vertices on 61 to those on 62 has to be at least three 
links. One can ensure this without going through a tedious check, by gluing cyclically 
together three copies of the original (spherical) manifold, i.e. 6i/6 2 , 6^/62, 6'//6 2 . The 
(double-) prime is used to distinguish the different copies. The boundaries 61, 6 2 are 
chosen such that they have no common vertex. 

The S 1 x S 1 x S 1 x S 1 manifold can be built out of the regular 3 4 square torus by 
dividing each elementary 4(i-cube into 4-simplices in the following way. For each cube 
we mark two points po = (0, 0, 0, 0) and p 4 = (1, 1, 1, 1) lying on the opposite ends of the 
main diagonal and connect them by one of the shortest paths going along the edges of the 
cube. The shortest path goes through 4 edges, each in a different direction, and through 
three points, say Pi,P2,P3- There are 24 such paths. Each set of points Po,Pi,P2,P3,P4 
forms a 4-simplex. 
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4 Finite Size Analysis and Results 



The elongated phase of simplicial gravity is well described as an ensemble of branched 
polymers §, §, §. 

This means that in the large volume limit one can use the Ansatz 

F(k 2 , iV 4 ) = NJ (k 2 ) + (7 - 3) logiV4 + /i(«2) (13) 

for the free energy. The corrections are expected to be of order 0(1/ N4). The correction 
coefficient 7 is assumed to depend only on the genus g of the underlying branched 
polymer structure. 

Differentiating ([13]) with respect to k 2 one sees that for the action density r (^) 

r = + ^ (14) 

the logarithmic corrections disappear. Therefore one should take the next corrections, 
namely I/N4, into account. They can appear in r for purely kinematic reasons. To un- 
derstand their origin consider the limit of large positive k 2 , in which only triangulations 
maximizing N 2 contribute to the sum Such triangulations can be obtained from 
barycentric subdivisions of 4-simplices applied successively to a minimal starting config- 
uration, possibly mixed with micro-canonical transformation, which do not change N 2 
and N4. By the minimal configuration we mean the minimal volume triangulation which 
maximizes N 2 . For the barycentric subdivisions one gets the relation N 2 = 5/2 N4 + c°, 
where the constant 

c ° = jV 2 ° - 5/2iV 4 ° (15) 

characterizes the initial minimal configuration (indicated by the index 0). The number 
N 2 of triangles is related to the action density r = N2/N4 = 5/2 + c° /N4. This means, 
the constant c° leads to I/./V4 corrections of the action density. The contributions to 
the sum (0) from triangulations built from non-minimal ones (ie smaller N 2 ) are sup- 
pressed exponentially by exp — n 2 (c° — c), where c characterizes the non-minimal start 
configuration with c(k 2 ) < c°. 

For the sphere, the minimal configuration is the surface of a 5-simplex, and therefore 
c° is known. For the other topologies we extracted the number listed below by a numerical 
experiment. With the standard topology preserving Monte Carlo process we use a cooling 
procedure, in which we increase k 4 , to decrease N4, and k 2 to maximize N 2 . In fact we 
have to increase k 2 slowly compared to K4, because increasing k 2 also increases the 
pseudo-critical coupling hi(n 2 ). 
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Figure 1: The volume dependence of the action density r for three different topologies 
in the elongated phase at k 2 = 2.0 



For the topologies studied, we found that the following configurations are minimal : 



(16) 





S 4 


N° A 


= 6, 


iV° 


= 20 


^c° = 


■5. 




x S 1 


iV 4 ° 


= 110, 


N$ 


= 44 


^c° = 


o. 


S 1 xS 1 x S 1 


xS 1 


N° 4 


= 704, 


N° 


= 1472 : 


^c° = 


28 



For the sphere S 14 the I/./V4 effect is very difficult to detect already for the volumes in 
the range of a few thousand 4-simplices and it would require extremely long runs to 
reduce the error bars below it. For manifolds S 3 x S 1 the effect is not present at all. The 
corrections are, however, two orders of magnitude larger for S 1 x S 1 x S 1 x S 1 and are 
measurable in the volume range used in the simulations. This estimation of the 1/N& 
effect is exact for infinite positive k 2 but one expects it to work although with a slowly 
varying coefficient c{k?) in the entire elongated phase. 

In figjl] we show the action density r measured in the elongated phase for n 2 = 2.0 
against I/N4. On the same figure we display the curve + c /A^ 4 with c° = 288, and 

= 2.482 which fits the data points very well. We note that that does not depend 
on topology, at least for those used in the simulation. 
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We find, with the statistics available, no volume- or topology-dependence of k 4 . This 
is compatible with the Ansatz 

« 4 («2,iV4)=/o(«2) + I T^, (17) 

because 7 is known to be 0(1). With the method used, the correction could only be 
separated from the statistical noise with an disproportionate amount of computer time. 
One could instead use the methods used in || to determine 7. We find = 5.659(4) 
for K2 = 2.0 in the infinite volume limit. 

In the crumpled phase, we use the power-law Ansatz 

F(k 2 ,N a ) =N i U{K 2 ) + h{K 2 )Nt (18) 

Taking the derivative with respect to N 4 one gets: 

R 4 ( K2} N 4 ) = / (k 2 ) + Sh^Nt 1 . (19) 

For the action density (|5|) one finds 

r = fifa) + fiMNt 1 (20) 
We checked the Ansatz by fitting our numerical data for r(N 4 ) (|20|): 



Topology 


s 


r°° = fi 


log(/{) 


xVdof 


S 4 
(S 1 ) 4 
S 1 x S 3 


0.5 (2) 
0.6 (2) 
0.5 (2) 


2.039 (+0.010 / -0.013) 
2.028 (+0.008 / -0.016) 
2.038 (+0.010 / -0.021) 


1.9 (1.3) 
1.0 (0.9) 
1.8 (1.3) 


0.78 
0.87 
0.31 



and for k a (N 4 ) to (0): 



Topology 
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*t = h 


log(5/i) 


X 2 /dof 


S 4 
(S 1 ) 4 
S 1 xS 3 


0.6 (2) 
0.6 (3) 
0.6 (2) 


1.20 (2) 
1.20 (2) 
1.20 (2) 


1.4 (1.3) 
1.4 (1.5) 
1.4 (1.3) 


0.22 
0.10 
0.18 



In the tables we give the logarithms of f[ and 8 fx, while these have approximately 
symmetric errors. One can see that the the values do not depend, within errors, on 
topology. In figure |2] we show the numerical data for K4 for the three topologies and 
8 = 0.5. 
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Figure 2: The volume dependence of the critical K4 for three different topologies in the 
crumpled phase at k 2 = 0.0 



The value 5 = 0.5 can be understood by looking at the typical configurations, that 
dominate the ensemble for large negative k 2 . These are configurations, which minimize 
the free energy, i.e. which have for a fixed volume the minimal number of vertices 
and thereby the minimal number N2 of triangles. For the three-dimensional case such 



configurations were constructed in []T0|| . This construction can easily be extended to the 
four- dimensional case. One starts with a 2c? triangulation of the sphere with t triangles. 
At each triangle one builds a 4<i-pancake neighborhood from q 4-simplices lying around 
the triangle in such a way that the links opposite to this triangle form a circle. An 
opposite link is defined to be a link, which does not share a common vertex with the 
triangle. The next step is to put neighboring pancakes together by identifying these 
circles with the 3d faces of the neighboring pancakes. Each pancake has three such 
faces which lie between the circle and an edge of the basic triangle. It also has three 
neighboring pancakes. After this step one gets an S A sphere with t * q 4-simplices and 
2 + 1 1 2 + q vertices. The highest connectivity N4 ~ (iVo — 2) 2 /2 is reached when t = 2q. 
The number N2 of triangles is a linear combination of Nq and N4. Using equation [| one 



finds r oc 1/ / 'sqrt(N^) ie. 5 = 0.5 for these configurations. We want to note, that in |TT 



an argument in favor of 5 = 0.75 was given. This can not be excluded by our numerical 
data. 



In [12] evidence was given, that for spherical topology the phase transition is of first 
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order. This was confirmed in [13| , |14] . We could observe flip-flops in the action density 
also for the two tori under investigation. We interpret this as a hint, that the transition 
is of first order in these cases as well. 

Finally we want to comment on the behavior of the algorithm in the crumpled phase. 
The typical configuration of simplicial gravity in this regime has one so-called singular 
link. The local volume of a singular link, i.e. the number of four-simplices which contain 
this link, diverges, when the volume of the entire configuration goes to infinity. The 
relaxation time, i.e. the number of Monte Carlo sweeps required until the singular link 
appeared in the configuration, was extremely long for all topologies. This effect is even 
more pronounced if one starts with a branched-polymer like configuration. For the (S 11 ) 4 - 
torus and small volumes (less than 64k 4-Simplices) the situation is even worse. It was, 
at least for the run-length used in our numerical experiments, impossible to relax to such 
singular configurations. On the other hand, we know that they exist, because they could 
be reached by shrinking down larger configurations containing a singular link. It seems 
it is difficult for the algorithm to deal with two different defects, the singular link and 
the hole, at the same time. We will discuss this point, which is directly related to the 
question of practical ergodicity, more carefully in a forthcoming publication. 



5 Discussion and Conclusions 

In this paper we have investigated the behavior of the entropy density and the curvature 
for three different topologies in four dimensional simplicial gravity. We employed man- 
ifolds consisting of between 4000 and 64000 simplices. We concentrated on two values 
of the gravitational coupling constant k 2 , one in the crumpled and one in the branched 
polymer phase. We found that in both the cases the value of the entropy density and 
curvature in the infinite volume limit are equal for these three topologies. This gives 
further support to the conjecture that these limits exist, a question which was discussed 
in the literature some time ago [jKJ, (TT|, [|TjJ. 



Furthermore, it gives support to the hypothesis that all topologies contribute to a 
sum over topologies, like in two dimensions. The value of the entropy exponent can be 



compared with the estimates given in |jl7| , based on a summation over all distributions 
of curvature. This estimate, which surprisingly enough is exact for the leading term in 
two dimensions, does, however, not directly give results in agreement with our numerical 
data for k 4 . 

We further analyzed in detail the finite size effects, and compared them to estimates 
from simple kinematic arguments. This approach explain very well the finite size cor- 
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rections. It may, of course still be possible that those effects are of a more complicated 
nature near the transition. 

Finally, as also seen in previous investigations, we observe that the algorithm has very 
long relaxation times in the crumpled phase. We even found that for the (Si) 4 torus, the 
ground state seemed not to accessible using the approximately fixed N± algorithm, but 
only by passing through states with much larger A/4- values. This breakdown of practical 
ergodicity is still under investigation. 
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